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NOMENCLATURE 


A nuclear mass number 

B(e) average slope parameter of nucleon-nucleon scattering amplitude, 

fm 

c Woods-Saxon surface diffuseness, fm 

e two-nucleon kinetic energy in their center of mass frame, GeV 

J defined in Eq. (21) 

m nucleon mass, kg 

R Woods-Saxon half-density radius, fm 

r position vector, fm 

r u radius of equivalent uniform distribution, fm 

t average two-nucleon transition amplitude, MeV 

t 0 defined in Eq. (3) 

V 0 defined in Eq. (16) 

v opt optical potential (defined in Eq. (15)), MeV 

x relative position vector of projectile, fm 

y two-nucleon relative position vector, fm 

z position vector of projectile in beam direction, fm 

°(e) average ratio of real part to imaginary part of nucleon-nucleon 
scattering amplitude 

3 defined in Eq. (6) 

6(r) Dirac delta function 

p nuclear density, fm -3 

P 0 normalization constant in Eq. (20), fm" 3 

o(e) average nucleon-nucleon total cross section, mb 



Subscripts : 

P projectile 

T target 

Arrows over symbols indicate vectors. 
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Lawrence W. Townsend 
Langley Research Center 

and 

Hari B. Bidasaria 
Old Dominion University 


SUMMARY 

A simple analytic method for calculating nucleon-nucleus optical 
potentials using a single-folding of a Gaussian two-body interaction with 
an arbitrary nuclear distribution is presented. When applied to 
proton-lead elastic scattering, the predicted real part of the Woods-Saxon 
potential is in substantial agreement with the experimentally-determined 
phenomenological potential, although there are no adjustable parameters. 
In addition, the volume integrals of both real potentials are nearly 
identical. 



INTRODUCTION 


In order to properly assess the optimum shield requirements for future 
manned space efforts, especially missions of long duration, accurate 
methods for predicting the interactions of cosmic rays with spacecraft 
materials and inhabitants are required. Although theories capable of 
accurately predicting high-energy nucleon-nucleus and nucleus-nucleus cross 
sections exist (refs. 1 through 6), disagreements between theory and 
experiment are evident at low energies. These are probably due, in large 
part, to the inapplicability of the eikonal approximation at low energies 
(refs. 1 and 2). In a recent work (ref. 7), substantial improvement in the 
agreement, at low energies, for light and medium nuclei was noted when 
partial wave analyses, utilizing complex WKB (Wentzel-Kramers-Brillouin) 
phase shifts, were performed. The method outlined in reference 7, however, 
appears to be suitable only for collisions involving light and medium 
nuclei whose matter densities are Gaussian or harmonic well distributions 
since their nuclear optical potentials can be analytically determined 
(ref. 8). For heavier nuclei whose density distributions are neither 
Gaussian or harmonic well shapes, their optical potentials are not readily 
reduced to analytic forms (ref. 1). In this work, a method is presented 
for obtaining an approximate analytic expression for the nucleon-nucleus 
optical potential, involving the single-folding of an arbitrary target 
density distribution with a Gaussian two-nucleon interaction. As such it 
represents an initial effort toward discovery of approximate analytic 
methods for eventual use in evaluating double-folding optical potentials 
involving arbitrary nuclear density distributions. When applied to 1 GeV 
proton-lead collisions, this method yields a predicted real part of the 
nuclear potential which is in excellent agreement with the available 
phenomenological results obtained from elastic scattering experiments 


2 


THEORETICAL ANALYSIS 


In previous work (ref. 2) it was shown that the optical potential 
approximation to the exact nucleus-nucleus multiple-scattering series is 


V(x) ApAj / d 3 z Pj(z) / d 3 y p p (x + y + z) t(e,y) 


( 1 ) 


(ref 6 ^ he . s two " bod y transition amplitude, averaged over constituent types 


t(e,y) = t Q exp[-y 2 /2B(e)] 


( 2 ) 


with 


t Q = -(e/m) ? a(e)[o(e) + i] [2irB(e)] 7 


(3) 


? nd p T , and P .P , are the target and projectile single-particle number 
n t h e . r f dens1t,es ); 1" equation (3), e is the two-nucleon kinetic energy 
n their center of mass frame, and a(e), a(e), and 8(e) are the usual 
nucleon-nucleon scattering parameters (refs. 1 and 2). 

For nucleon-nucleus scattering involvinq a taraet nurlpn<; nf macc 
number A T , the projectile single-particle density is 


P p (x + y + z) = 6(x + y + z) 


(4) 


with A P = 1. Equation (1) then reduces, for nucleon-nucleus scattering. 


^opt^ x ) f d ^ z t(e, x + z) 

Inserting equation (2) into (5) and letting 
6 2 = [ 28 ( e )]- 1 
yields 

V opt (x ) = 27rt 0 A T ^ z2(Jz p t(z) exp[-B 2 (x 2 + z 2 )] 
0 

/° exp[-2xz e 2 cos 0 ] sin 0dG 

77 


(5) 


( 6 ) 


(7) 
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The angular integration is of the form 


/"exp(-q cos ©) sin Ode = [exp(q) - exp(-q)]/q (8) 


which, upon collecting exponents, yields 

V op t(*> = ( 7lt 0 A T /e2 x ) /" zdz p T ( z ) <exp[-3 2 (z - x) 2 ] 

0 

-exp[-B 2 (z + x) 2 ]} 

Evaluating the integrals in equation (9) gives 

/"zdz p T (z)exp[-3 2 (z - x) 2 ] = 

0 

,-l 


and 


f / (x + 3’ 1 s)p t (x + 3"^ s)exp(-s 2 )d; 
-3x 1 


/"zdz p T (z)exp[-B 2 (x + z) 2 ] = 


3 * P (3 * s - x) p t ( 3 - ^ s - x)exp(-s 2 )ds 
Sx 


(9) 


( 10 ) 


( 11 ) 


Incorporating the results from equations (10) and (11) into (9), and 
assuming that the density distribution is spherically symmetric yields, 
after some algebra, 

V opt^ = irt o A T 0 " 3 { ( 6 x ) _1 fsds exp(-s 2 ) [p T (x + 3 _1 s) 

- p t (x - 3” 1 s)] (12) 

+ /" ds exp(-s 2 ) p t (x + 3 _1 s)} 


The integrals in equation (12) are of the same form as those in 
equation (86) of reference 1, where a similar method was utilized to 
extract matter densities from nuclear charge densities. 
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The first integral, is generally smaller than the second since it 
contributes only when x is near the nuclear edge. For the extreme case 
where pj is a finite uniform distribution, the ratio of the first 
integral to the second, at the uniform radius, r u , is 


Ratio = r _1 [2B(e)/ir] 


1 


(13) 


Noting that this ratio is a maximum for large B(e) and small r u , which is 
the situation for light target nuclei and high energies, an estimate of the 
ratio is made. Choosing r u = 3 fm (a nominal value for a lithium nucleus 
from reference 10) and B = .5 fm 2 (ref. 11) gives a maximum ratio of 19 
percent. It is expected that the error resulting from neglect of the first 

integral, for realistic densities, will be substantially less than this 
v & 1 ue . 


Neglecting the first term in equation (12) yields an approximate 
nuclear optical potential 


V opt( x ) " 7rt 0 A T B " 3 ^ ds ex P(-s 2 ) p t (x + 6' 1 s) 


(14) 


For an arbitrary density distribution, the integral can be approximated by 
a two-point Gauss-Hermite quadrature formula (ref. 12) to yield an analytic 
nuclear optical potential 3 

V opt( x ) = ( v 0 A T /2 Hp t (x + /B) + p T ( x - /B)] (15) 


with 


1 

v 0 = -(e/m)^ a(e)[a(e) + i] 


(16) 


Should the need arise to include the first integral in equation (12) it 
can be similarly approximated by changing variables such that 

/” sds exp(-s 2 )[p T (x + ef 1 s) - p T ( x - e' 1 s)] 


= j /“dp exp(-p)[p T (x + (p/e 2 ) 7 ) ( 17 ) 


_ 1 

- p t (x - ( p/6 2 ) 7 )] 


formula 1 (ref th 12) nte9ral t0 ^ eaSl1y a PP roximat ed by a Laguerre quadrature 
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RESULTS FOR PROTON-LEAD ELASTIC SCATTERING 


To illustrate the application of equation (15), the real part of the 
nuclear optical potential for proton-lead elastic scattering at an incident 
kinetic energy of 1.04 GeV is calculated. Thus, taking the real part of 
equation (15) gives 


1 

V (?) = a < e ) «( e ) V 2 Hp T ( x + ^ B ) + p y( x - V*B)] (18) 


which, for 1.04 GeV protons colliding with lead yields 


V opt (x) = 9641.5 [p t (x + .49) + p T (x - . 49 )] 


(19) 


For pj given in fm~ (normalized to unity), 
given in MeV. Choosing py of lead to be a 
using the methods in references 1 and 2 yields. 


Vopt( x ) 1n equation 
Woods-Saxon charge 


(19) is 
density 


p y(r) = P Q tl + exp [ (r - R)/c]} -1 


( 20 ) 


where R = 6.624 fm, c = 0.860 fm, and p 0 = 7.059 x lO" 4 fm -3 . The 

results obtained from equations (19) and (20) are plotted in figure 1. 

Also plotted in figure 1 are the phenomenological Woods-Saxon potential 
results obtained from elastic scattering experiments (ref. 9). The 

agreement between theory and experiment is impressive particularly since 

there are no arbitrarily adjustable parameters in the theory. In addition, 
the volume integrals, J , of both real potentials 

J = 4tt J°° r 2 V(r) dr (21) 


agree te within 0.02 percent. 

Finally, we note from equation (13) that the maximum expected error 
due to the neglect of the first integral in equation (12) is less than 6 
percent. The actual error, for the density used in the calculations 
(eq. (20)), was found by numerical methods to be less than 5 percent 
(typically 2-3 percent) for values of x up to 12 fm. 



CONCLUDING REMARKS 


In this work a simple analytic method for approximating nuclear 
optical potential integrals involving the single-folding for a Gaussian 
interaction with arbitrary nuclear distributions was presented. Applying 
the method to 1.04 GeV proton-lead collisions, the real part of the 
predicted potential was found to be in remarkably good agreement with the 
phenomenological results obtained from elastic scattering experiments. The 
resulting Woods-Saxon radial shape, for the potential, obtained when a 
Woods-Saxon target density was utilized also suggests that there may be 
some credibility in the usual assumption that the spatial dependence of the 
nuclear potential should follow that of the associated nuclear density 
distribution. It is anticipated that the analytic methods described herein 
will be useful for determining complex WKB solutions to low energy 
scattering problems and for other nuclear potential calculations where 
purely numerical methods may not be appropriate. 
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Figure 1. - Real part of the proton - lead optical potential at 1.04 GeV incident 
kinetic energy. 
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